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We have extended through the high-temperature expansion of the second field derivative 
of the susceptibility for Ising models of general spin, with nearest-neighbor interactions, on the 
simple cubic and the body-centered cubic lattices. Moreover the expansions for the nearest-neighbor 
correlation function, the susceptibility and the second correlation moment have been extended up to 
. Taking advantage of these new data, we can improve the accuracy of direct estimates of critical 
exponents and of hyper-universal combinations of critical amplitudes such as the renormalized four- 
point coupling Qr or the quantity usually denoted by R^. In particular, we obtain 7 = 1.2371(1), 
V = 0.6299(2), 74 = 4.3647(20), gr = 1.404(3) and = 0.2668(5). We have used a variety of series 
extrapolation procedures and, in some of the analyses, we have assumed that the leading correction- 
to-scaling exponent 9 is universal and roughly known. We have also verified, to high precision, the 
validity of the hyperscaling relation and of the universality property both with regard to the lattice 
structure and to the value of the spin. 

PACS numbers: 05.50.-|-q, 11.15.Ha, 64.60.Cn, 75.10.Hk 


I. INTRODUCTION 

The numerical study of the critical properties of the spin-S' Ising models with nearest-neighbor interactions hacLap. 
important historical role in the chain of arguments leading to the modern formulation of the universality hypothesisEfu 
for the critical phenomena and, in particular, to the concept of universality class. It was in a study of the susceptibility 
x(/3; S) for the general-spin Ising models omthe face-centered-cubic (fee) lattice by high-temperature (HT) expansions 
through /3®, that C. Domb and M.F. SykesQ first pointed out that the exponent 7 , which characterizes the divergence 
of Xj was rmadaly independent of S and guessed for it a universal “daltonian” value 7 = 5/4. Later on, when 
longer seriesclLj both for the fee and for other lattices were derived, a weak dependence of the exponents 7 and iz 
on S emerged from the HT analyses, but was soon correctly ascribed to the occurrence of non-analytic “confluent 
corrections to scaling” (CCS) rather than to a failure of the universality. Inj-thp^ years the existence of CCS had|t>epi| 
inferred by various authors both from the numerical analysis of HT serieidril3 and from phenomenological htsEJTH 
to high precision experimental data for some systems close to criticality. Eventually the status of the CCS was more 
firmly establishedliZi in the context of the Renormalization Group (RG) theoryEl. It was therefore recognized very early 
that the accurate determination of the critical exponents in numerical or experimental studies and, as a consequence, 
the feasibility of stringent verifications both of the universality hypothesis and of the scaling and hyperscaling relations 
require a close control over the GGS. For many years, however, in two and in three dimensions, HT series were available 
only for few observables and, generally, they were barely sufficient to conject.«^|-tlpi presence of GGS, but definitely 
too short to make a numerically accurate discussion of these models possibleSMEJlij. Still presently, the expansions 
of xiP'j S) and of the second moment of the corijelation function pL 2 {P; S) on the simple-cubic (sc) lattice, for spin 
S > l/2_can be found explicitly in the literatures only up to the order (3^'^. The data files by R. Roskies and P.M. 
Sackettc3 made an extension of these series through feasible for the sc and the bcc lattices, rbut did not drastically 
changeJffie situation. On the fee lattice, the HT series initially derived through in Ref.0, were later extended 
in Ref.EJ to order (3^*. Fortunately, in the case of the body-oentered-cubic (bcc) lattice, decisive progress occurred 
already two decades ago, with the computation by B.G.NickeEd of exmnsions for x(/3; S) and fj.2if3; S) through /3^^. 
(To our knowledge, only the series for S = 1/2,1, 2, 00 were published^.) r—■ r—. 

By allowing to some extent for the leading GGS, the first modern analyses of the extended bcc seriestdcS improved 
significantly the accuracy in the verification of universality with respect to the magnitude of the spin. Moreover, in 
the mentioned studies (as well as in later analysesEj mainly devoted to the S' = 1/2 case), the central estimates of the 
susceptibility and the correlation-length exponents were reduced up-ta-« 1% with respect to the values 7 = 1.250(3) 
and v = 0.638(2), initially guessed in RefE and later confirmedEilS by various studies. This development also 
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contributed to.yttleH a long-standing controversy raised by the results of Refsll^’tSTEj, which stimulated the studies 
of Refs.L3~EjHE3, on the validity of hyperscalii^a^, more generally, on the consistenc;^f the results from the HT 
analyses with the ccpje^onding RG estimateaS’cJ^, either in the e-expansion approach^ or in the hxed-dimension 
perturbative schemeCjO. 

One should also note that for the second held-derivative of the susceptibility X4(/3; S) and for the nearest-neighbor 
correlation function G(/3; 5'), the published data are even less abundant. On the sc lattke, series for X4(/3; S) can be 
derived from the data files of RefEu up to cpier and up to from the data of Ref.o fordiie bcc lattice. On the 
fee lattice, series for X4(/3;S') are available^ through (3^^. For general spin, only expansionsaO of G{P;S) through 
on the fee lattice have been published. A summary of the HT expansions available until now for the Ising models 
of general spin appears in Table |. 

We have ^en pursuing a long-term pspiact to improve the algorithms and the codes for HT expansions in two- 
dimensionalEil and in three-dimensionalEuca lattice spin models, keeping up with the steady increase of computer 
performances, and periodically updating the numerical analjjppwhenever we could significantly extend the series. By 
using an appropriately renormalized linked-cluster methodaruM, we have now added from four up to thirteen terms 
to the HT expansions for various observables of the general spin-S* Ising models on the sc and the bcc lattices. In 
this paper we shall examine the expansions of x(/3; S) and /r2(/3; S) up to order and of X4(/3; S) up to on both, 
lattices. These data have been derived by slightly improving the thoroughly tested code which recently producedE£l 
our series through for x(/3; 1/2) and 1/2) on both lattices. The extension of the series is by far the hardest 
part of this work, but we will not enter enter here into the details of our procedure. To give an idea of the required 
computational effort, it will suffice to mention that our improved codes take minutes of CPU time on a COMPAQ 
Alpha XPIOOO (500 MHz) single-processor workstation to reproduce the known series through whereas several 
days are necessary to add the following four orders. From the graph-theoretical point of view, it is the expansion of X 4 
throughA?^^ which involves the most laborious part of the calculation: in the simplest vertex renormalized expansion 
schemeEJ it would require the generation and the evaluation of over 10® topologically inequivalent graphs. However, 
devising a careful strategy of in-depth renormalizations, the expected size of the calculation has been reduced by at 
least two orders of magnitude. On the other hand, from a purely computational standpoint, the calculation of the sc 
lattice constants for the second moment of the correlation function is the most demanding part of the job in terms of 
CPU-time. 

The correctness of our codes is ensured by numerous internal consistency checks, as well as by their ability to 
reproduce established results already available in simpler particular cases, such as the square-lattice two-dimensional 
spin 1/2 Ising model or the one-dimensional spin-S' Ising models. Of course, our codes also reproduce the old 
computation of Ref.ufor S = 1, 2, oo on the square and the bcc lattices, and, as far as there is overlap, also the recent 
computation of Ref.H for S' = 1/2 on the bcc lattice. 

Using this vast library of partially new high-order series data and in particular our significantly extended series for 
X4(/3f-S)Tnwe can resume from a vantage point the very accurate studies performed on series 0(/3®^) for x and /i 2 in 
RefsE3nH3 and present an even more extensive and detailed survey of tha critical behavior for the spin-S Ising models. 
In spite of the remarkable advances achieved by the calculations of Ref.E3cZl which removed away from the foreground 
the universality and the hyperscaling issues, further extensions of the HT data still remain of great interest. They are 
instrumental in the continuing efforts to gain a higher accuracy in the estimates of the critical parameters and, more 
generally, to perform more stringent tests of hyperscaling and of universality, with respect both to the value of the 
spin and to the lattice structure. These are certainly welcome results, since it is fair to say that the actual verification 
of such basic properties is still only moderately accurate, although no doubts persist anymore about their validity. 
Of course, one must be aware that the computational complexity of the calculation of higher-order series coefficients 
grows much faster than the precision in the evaluation of the critical parameters that can be obtained from them by 
the presently available numerical tools. Therefore the higher-order computations should be accompanied also by an 
effort to improve the techniques of analysis or, at least, by a careful comparison of the results obtained by a variety 
of methods. 

The paper is organized as follows. In the next Section, we set our notations and definitions. In Sec.HI we state 
the assumptions underlying our axialysis and its aims. The numerical procedut^ we have used, namely the modified- 
ratio methods introduced in Ref.O or the differential approximant methodsSEu, as well as the corresponding results 
of the series analysis, are discussed in Sec. IV-VHI. In Sec. IX we compare our estimates with those of the most 
recent literature. The last few Sections present our results for the critical amplitudes of the observables that have been 
expanded and for some (hyper)-universal combinationsEil of these amplitudes. In order to make our analysis completely 
reproducible and to provide a epavenient source of data for further work, without overburdening this paper, we have 
collected into a separate reported, available on request, the complete expansions of the nearest-neighbor correlation 
function, of the susceptibility, of the second moment of the correlation function and of the second held derivative of 
the susceptibility for spin S' = 1/2,1, 3/2, 2, 5/2, 3, 7/2,4, 5, oo, on the sq, sc and the bcc lattices. 
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II. THE SPIN-S' ISING MODELS 


The spin-S' Ising models are defined by the Hamiltonian: 

H{s} = ^ s(x)s(f') (1) 

where J is the exchange coupling, and s{x) = s^{x)/S with s^{x) a classical spin variable at the lattice site x, 
taking the 25” + 1 values —S', —S + 1,..., S — I, S. The sum runs over all nearest-neighbor pairs of sites. We shall 
consider expansions in the usual HT variable (3 = J/ksT called “inverse temperature” for brevity. 

In the high-temperature phase, the basic observables are the connected 2n-spin correlation functions. Here we shall 
limit our study to quantities related to the two-spin correlation functions {s{x)s{y))c and to the four-spin correlation 
functions {s{x)s{y)s{z)s{t))c- 

In particular, we shall consider the nearest-neighbor correlation function 

OO 

G{p-,S) = {s{M5))c = Y.hriSW ( 2 ) 

r=0 


where d is a nearest-neighbor lattice vector. 

The internal energy per spin is defined in terms of G'(/3; S) by 

U{(3;S) = -^G{(3;S) 

where q is the lattice coordination number. 

The specific heat is the temperature-derivative of the internal energy at fixed zero external field 


GH{P;S)/kB 


ql3‘^ dG{f3; S) 
~2 


In terms of x(/3; S), the zero-field reduced susceptibility, 


OO 

x(/3; S) = ^(s(0)s(f))c = Cr{S)(3^ 

X r—0 


and of /r2(/3; S), the second moment of the correlation function. 


( 3 ) 


( 4 ) 


( 5 ) 


d2{(3\S) = ^x2(s(0)s(x))c = Y‘^r{S)f3^ 


the “second-moment correlation length” ^{(3;S) is defined by 

t^2{P;S) 


eiP;S) = 


= YtriSW- 


W-,S) 

The second field-derivative of the susceptibility X4(/3; S) is defined by 

OO 

X4(/3; S)=Y (s(0)'S(a;)s(2/)s(2:))c = Y ^r{S)l3G 

x,y,z r—0 

Notice that these definitions ensure the existence of a non-trivial limit as S' ^ oo. 


( 6 ) 


( 7 ) 


( 8 ) 
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III. ASSUMPTIONS AND AIMS OF THE SERIES ANALYSIS 


In the universality class of the spin-S' Ising models, the asymptotic behavior of the susceptibility as /3 ^ Pti^) 
from below, is expected to be 

X* (/3; S) ~ C* {S)t* {l + a* iS)T* {Sf + ... + b* (5)t# (5) + ...) (9) 

where r'^iS) = 1 — (3/(if {S) is the reduced inverse temperature. We have introduced here the superscript # which 
stands for either ^ or bcc, as appropriate, and will be used hereafter only when useful. Eq.(H), often called the 
Wegner expansiorO, specifies how the dominant scaling behavior, characterized by the universal critical exponent 7 
and by the critical amplitude C'^(S'), is modified by analytic and nonanalytic confluent corrections to scaling(CCS) in 
a close vicinity of the critical point. The leading nonanalytic CCS is characterized by a universal exponent 9 and by 
an amplitude a^(S'). The critical amplitudes (S), (S), (S), as well as the inverse critical temperature f3f{S), 

are nonuniversal: namely they depend on the spin S and on the lattice structure, as stressed by the notation. The 
analogous asymptotic behaviors of the correlation length 

ew; s) f*iS)r*iS)-''(l + afiS)T*iSf + ... + bfiS)r*iS) + ...) (10) 

of the specific heat 

C*{(3; S)/kB A*{S)t*{S)-^{i + a*(S)T*{Sf + ... + b*{S)T*{S) + ...) (11) 

and of X 4 (/ 3 ; S) 

xf (/3; S) ~ -Ct{S)T*{S)-^^ (1 + af (5)r#(5)^ + ... + bt{S)T*{S) + ... ) ( 12 ) 


as well as of the other singular observables, are characterized by different critical exponents and by different (nonuni¬ 
versal) critical amplitudes etc., but all contain the same leading confluent exponent 6. Notice that we 

have freely chosen in equations d)-® between the conventions of Ref.l^ and those of Ref.l^, since the notation for 
the amplitudes is not yet completely standardized. 

Usually the exponent 74 is expressed in terms of 7 and of A, the “gap” exponent associated with the critical 
behavior of the higher field-derivatives of the free energy, as follows: 74 = 7 -|- 2A. Here A = /3 -|- 7 and (3 denotes 
the magnetization .eTOonent only in this formula and in the scaling relation to be quoted before eq.(|20|). From RG 
calculationscJoa’y it is expected that 0 ~ 0.5 for the Ising universality class. 

For later use, we observe that, if the singularity closest to the origin in the complex (3 plane is the critical singularity, 
then eq.(^) implies the following asymptotic behavior for the expansion coefficients cf{S) of x^(/3; S) 




c#iS) = C*iS)^p--{S) 


1 -f 


r(7) «f(^) 

r(7 — 6 ) 


0 (l/n) 


(13) 


In the Ising universality class, the expected 0(l/n^®) contributions in eq.(p^ should be practically degenerate with 
the analytic 1/n corrections. For bipartite lattices- the higher-order corrections in eq. include terms 0 (l/n^+^ “) 
with alternating signs, which reflect the presenceEALS of a weak “antiferromagnetic” singularity at /3 = —Pf{S) with 
exponent 1 — a. Analogous formulae for the asymptotic behavior with respect to the order can be written for the 
expansion coefficients t^{S) of ^^(/3; 5)^ and ejf (5) of xt (z®! S). 

In terms of x? ^ and X 4 ) a “hyper-universal” combination of critical amplitudes denoted by gr and usually called 
the “dimensionless renormalized coupling constant”, can be defined in d dimensions by the limiting value of the ratio 


g*{fi;S) = 


3v*xt{S-,P) 

\Qtt^*{S]I3Yx*{S\I3Y 


(14) 


as t#( 5') ^ 0-1-. Here denotes the volume per lattice site (in 3 dimensions = 1 and =|A/3-\/3) and the 
normalization factor is chosen in order to match the conventional field theoretic definition of g^. We shall call 
g"^{f3-,S) the “effective renormalized coupling constant” at the inverse temperature (3. 

By eqs. (D-®: g'^{S;P) behaves as 
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g* (/3; S) ~ grT* (5)« + ...) 


(15) 


when t'^(S') ^ 0 +, with 


16 Trf#{SYC*{S)^ 


The Gunton- 


Buckinghamii~Ei 


inequality 


(16) 


'y + diy — 2A > 0 


(17) 


together with the Lebowit^l inequality xtiP', S) < 0, ensures that g'^{P; S) remains bounded and non-negative as 
T^(S') ^ 0-I-. The vanishing of g'^iPf — 0-I-; S) is a sufficient condition for Gaussian behavior at criticality, namely for 
the vanishing of the four-spin and oirthe higher-order connected correlation functions. In lattice field theory language, 
this corresponds to the “triviality”BlZa, numerically observed when d ^-^nd proved when d > 4, for the continuum 
field theory defined by the lattice model (with ferromagnetic couplingsQO) in the critical limit. If the inequality & 
holds as an equality 

X + di/ — 2A = 0 (18) 

(called the ‘hyperscaling relation’), if there are no logarithmic corrections to the scaling behavioi0E§ and if xt 
is nonvanishing, we have: 


g*{p-, S) c^gr{l + a*{S)r*iSf + ... ) (19) 

namely the “ effective coupling” g^(/3; S) tends to a universal nonzero limiting value gr as t#(5') 

Using the Essam-FisherLETEj scaling relation a -|- 2/3 -|- 7 = 2, eq.(|T^) can be rewritten as a relatiorO’o between a 
and V 


a = 2 — dv. ( 20 ) 

As a consequence, also the following combination of the critical amplitudes A^{S) and 

Rt = ( 21 ) 

is hyper-universal, as pointed outB by Stauffer, Ferer and Wortis. . 

The ratios a^/a^ pid^ 4 /ax) etc., of the amplitudes of the leading GGS are less studied!^, but not less 

interesting, universalEATJ critical observables. 

In the rest of this paper we shall employ our HT series to estimate the critical parameters defined by eqs. i)-®, 
(H) and (| 2 ^), to check the validity of eqs.([l 8 |) and (| 2 ^) and of the universality property with respect to the value of 
the spin and to the lattice structure. 

In th^aetpal numerical analysis of finite-order HT expansions, the presence of the CCS will generally become 
manifest0c3c3 by small apparent violations both of the hyperscaling relations and of the universality properties, 
namely by a weak apparent dependence of the universal quantities on the lattice structure and on the value S of 
the spin. We shall point out this fact by explicitly indicating, in our numerical estimates of the universal quantities, 
the value of the spin and the lattice structure for the series used in the analysis. For instance, will denote 

the numerical estimate of the universal exponent 7 obtained from the series x^‘’‘’(/3; S'). This notation will help to 
emphasize how small the mentioned effects of apparent nonuniversality are reduced if the HT expansions can be 
pushed to a sufficiently high order, provided that the numerical tools of the analysis can, at least approximately, allow 
for the leading CCS. 

Part of our analysis will rely upon the main assumption that the exponent d of the leading CCS is universal and 
roughly known. A recent accurate RG recalculation of universal critical dataMI^ predicts the value 9 = 0.504(8) in the 
fixed-dimension perturbative approach, while within the e-expansion scheme, the updated estimate is 0 = 0.512(13). 
In the rest of this paper, we shall adopt as a reference value the fixed-dimension RG estimate 9^^^ = 0.504(8), when 
computing the central values of critical parameters by procedures biased with 9. Even if one has no compelling reason 
to supj^e that the uncertainty of the RG prediction of 9 is largely underestimated, (but this possibility is advocated 
in Ref.E3), the reliability of the 0-biased analyses presented here will be greater whenever their results are not too 
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sensitive to the precise value of 0. In the following Sections, it will be clear that, in most cases, we can tolerate an 
uncertainty of this exponent even several times larger than above indicated. We will show that most of our estirpate 
biased with =p0. 504(8) will be compatible also-adth h^er values, such as 0 = 0.52(3), proposed in Refs.^^, 
0 — 0.53(2) in Ref.N or even 0 = 0.54(3) from RefsE3’LJO~H. We can add that, both our direct HT evaluation of 0 
and the part of our series analysis which is not biased with the value of the exponent 0 , will be completely consistent 
with the above assumption. 

If appropriate, we shall provide detailed information on the 0-biased numerical results reported in our tables for a 
given quantity P, by indicating together with the central estimate, also the derivative dP/d0 evaluated at the reference 
value chosen for 0. Similarly, in the cases where the parameter estimates are biased with the value of a critical inverse 
temperature and/or of a critical exponent, for instance 7 , we shall report the corresponding derivatives 

dP/dj3c and/or dP/d^ computed at the specified reference values. As an example, for the critical amplitude of the 
susceptibility C'^(S'), our final estimate can be read as 

C*{S){erTor) + (9C#(5)/9/3,)(/3, - /S^^) + {dC*{S)/- 7 ™^). (22) 

Here both the estimate and its derivatives are evaluated for sharp values of j3c = and of 7 = 7 ’’®-^ and the 
error attached to the hrst term does not allow for the uncertainty of the bias parameters. Since the above expression 
describes how the central estimate of C^{S) changes under small variations of the bias parameters, comparisons with 
previous results in the literature, often based on slightly different assumptions, are made straightforward. 

As a final general remark, it is worth to mention that, due to the higher coordination number of the lattice, the bcc 
series approach their asymptotic structure eq. (H) generally faster than the sc series. For this reason, the bcc series are 
usually observed to yield more accurate estimates of the critical parameters than the corresponding sc series with the 
same number of coefficients and are often said to have a greater “effective length”. This fact will be confirmed here, 
and will be one of the reasons to draw our final best estimates from the analysis of the bcc series. Nevertheless, the 
sc lattice series remain very interesting, in particular because the non-universal informations obtained from them are 
directly comparable to the data from simulation studies, traditionally performed on the sc lattice. It is also interesting 
to notice that, for both lattices, the asymptotic behaviors set in more slowly in the most widely studied S = 1/2 
case. This is not surprising since the number of degrees freedom per site is proportional to the magnitude of the spin. 
A slower convergence is observed also for the higher moments of the correlation function, since in their construction 
larger weights are given to the correlations between farther sites for which the expansions are effectively shorter. As 
a consequence, on both the sc and the bcc lattices, the expansions of ^ shows a slower convergence than those of %. 

Basing on the assumptions above indicated, our analysis will aim to exhibit, within the family of the spin-S” 
Ising models, some consequences of the universality property, of the scaling and of the hyperscaling laws for the 
correlation functions, for the exponents and for various universal combinations of critical amplitudes. In particular 
our accurate verificatLon of the universality property will strengthen the justification xif a technique advocated long 
ago by J. Zinn-JustinCJ and independently by J.H.Chen, M.E. Fisher and B.G.NickeEj to improve the precision in 
the computation of the universal critical parameters of the Ising model. These authors argued that numerical study 
should address appropriate families of spin models parametrized by a continuous auxiliary variable and belonging to 
the same universality class as the Ising model. For specific values of this variable it is possible to select representative 
models for which the amplitudes of the leading confluent corrections to scaling are negligible and, as a consequence, 
the determination of the universal critical parameters can be more accurate. By relying on a similar prescription, we 
will also obtain very accurate estimates of some universal critical parameters. 


IV. ESTIMATES OF THE CRITICAL POINTS 


In this Section we shall examine the HT expansion of the susceptibility in zero field, for several values of the spin S 
on the sc and the bcc lattices. The series coefficients of the susceptibility generally show a very smooth dependence 
on the order of expansion and a relatively fast approach to their asymptotic forms. Therefore they are best suited to 
an accurate determination of the critical temperatures. The estimates so obtained will also be adopted to bias the 
calculation of the critical exponents and amplitudes. 

As we have already argued-in our previous studyl^ of the HT series for the S' = 1/2 case, the modified-ratio method 
introduced by J. Zinn-JustinEJ, (see alsoO), can lead to estimates of the critical inverse temperatures with an accuracy 
comparable or sometimes higher than the traditional differential approximant (DA) methods. Perhaps, the potential 
of this tool has not been properly appreciated, because so far it could not be used with series long enough. 

The method consists in evaluating (3f{S) from the approximant sequence 


(/3f (^))n = 


^n—1 


Sn + Sn-2 j 

2>Sn{Sn Sn— 2 ) 


(23) 
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with 


= 


^n-2 
CnCn—4 


) ^ -\-ln{ 


^n-3 


Cn—1^71—5 


)-i)/2. 


(24) 


Since the expected value of 9 is very nearly 1/2, by using eq. dil), the asymptotic behavior of the approximant 
sequence can be expressed as follows: 


(/3f(5))„=/3f(5)(l 


2r(7-0) n^+0 


(25) 


It is interesting to observe also that, if 0 = 1/2, the coefficient of the 0(l/n^) correction in eq.(^) is equal to 
(0^(5"))^ times a very small positive factor. Moreover, the coefficient of the 1/n-term in eq. dil) enters into eq.(|^) 
only at next higher orders. Since 0^(5") is expected to be small (though generally not negligible), these remarks help 
to understand how this method works and why it is much more efficient than the conventional ratio prescriptions. We 
can say that eq.(^^ provides an estimate of the leading “finite-order effects ”, namely of the corrections due to using 
series of finite length n. These are strictly analogous to the well known ’’finite-size effects” which have to be carefully 
considered to improve the data from simulations of finite systems. At the orders of expansion presently available 
eq.(|^) has already a reasonably accurate quantitative meaning. 

Although devised specifically to deal with the expected structure of the singularities, the procedure we have sketched 
is unbiased: namely no additional accurate information on other critical parameters must be used together with the 
series in order to get the estimate sequence. However, at the present orders of expansion, the n—dependence of {PciS))n 
is not saturated and, for sufficiently large the successive estimates show an evident residual trend, very nearly linear 
on a l/n^+® plot, as expected from eq.(E^. Small odd-even oscillations are superimposed to the main monotonic 
trend as a consequence of the above mentioned antiferromagnetic singularity (see the comments to eq. (HI))- These 
observations suggest that one can do something better than taking the highest-order available term of the sequence 
eq.(|2|) as the final estimate of Pf{S). The most obvious improvement consists in using the assumed known value of 9 
to fit the asymptotic behavior of the sequence and in taking the extrapolated value of the sequence as a better estimate 
of (3f{S). As usual, one should separately extrapolate to large n the odd and the even subsequences of {Pc{S))n, in 
order to deal properly also with the oscillations due to the antiferromagnetic singularity. Our extrapolation will be 
based on the successive pairs of terms in the approximant subsequences. Eventually a further minor adjustment of the 
results might be performed by a second (purely visual) extrapolation in order to allow also for a very small residual 
curvature of the plots due to the higher corrections in eq.(^). For instance, in the S' = 1/2 case on the sc lattice, the 
highest-order estimate from the extrapolation of the last pair in the odd-approximant subsequence is 0.22165646. In 
order to allow for the small residual curvature of the extrapolation sequence, this figure should probably be slightly 
reduced, to yield our final (and very conservative) estimate /3/^(I/2) = 0.221655(2). 

The set of our estimates for (if{S) is reported in Table ||. The errors we have indicated are small multiples (2-4) 
of the differences between the extrapolations of the two highest-order pairs of terms in the odd subsequences. In the 
same Table, we have also reported d(if{S)/d9 evaluated at 0 = 0”®-!'. As shown by our data, the above mentioned 
uncertainty in the value of 0’'®-^ turns out to be unimportant in the whole procedure, because it contributes only a 
small fraction of the final uncertainty of the estimates. 

In order to give an idea of the qualitative features of the method, for each value S of the spin examined in this 
study, we have plotted in Fig. |^the corresponding “normalized” approximant sequence (/3/°'’(5'))„/A^(5') vs. l/n^+®. 
We have taken the average of the extrapolated values of the even and the odd subsequences as the normalization 
factor N{S), introduced only to make the various plots easily comparable and conveniently fit all of them into a single 
figure. We have drawn as continuous lines the extrapolants of the last odd pair of terms in the sequences, whereas 
the dashed lines indicate the extrapolants of the last even pair. The difference between the extrapolated values of the 
odd and the even normalized subsequences, which is generally very small (for instance it is ~ 10“® in the bcc lattice 
case and at most four times as large in the sc case), provides a first rough indication that the oscillating corrections 
due to the antiferromagnetic singularity give only a small contribution to the uncertainty of the results. The final 
relative errors reported in Table y are generally much larger. 

If we refer to eq.(|^), the plots in Fig. strongly suggest that a^“(l/2) < 0, whereas > 0 for 5” > 2. Since 

|a^'’’'’(l)| and |a^'’’'’(3/2)| are very small, we cannot yet be completely sure about their sign. The smallness of these 
confluent corrections is confirmed observing that |9/9/°'’’(l)/90| and |9/3/°'’’(3/2)/90| are much smaller than for the 
other values of S. 

Simple model series with a structure specified by eq.(p^ (including the antiferromagnetic oscillating corrections) 
can mimic rather accurately the behavior of the spin-A Ising series for sufficiently high orders. Therefore numerical 
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experimentation with these model series can give us some intuition on the virtues and the limitations of the modified- 
ratio method and help to assess its accuracy. These tests add further confidence on our estimates of the relative error 
of Pc- On the other hand, it may take series significantly longer than those presently available to determine a^(S') with 
a precision better than a few percents, since the slopes of the approximant sequences provide only “effective” values 
of these amplitudes due to the residual influence of the higher-order corrections. Actually, the relative uncertainty of 
a^(S') can be larger, particularly so if its absolute value is very small. 

The sc lattice series have been studied in the same fashion and the results are illustrated by Fig. The main 
difference with respect to the bcc case is that all approximant sequences are decreasing, so that a^{S) < 0 for all S. 
It is also clear that, for this lattice, the rate of convergence of the approximant sequences to their asymptotic behavior 
is distinctly slower than in the bcc case. 

In order to gain further confidence in the estimates by the modified-ratio method, we must confirm at least their main 
features also by numerical tests of a different nature or involving different assumptions, thus reducing the probability of 
being misled by only apparent convergence. We have therefore performed also a more traditional unbiased analysis by 
first and second-order inhomogeneous DA’s yielding values of the critical inverse temperatures in essential agreement, 
to within their uncertainties, with those obtained from the modified ratio-method. In the case of the bcc lattice, 
for spin S = 1/2, the highest-order available DA estimates are slightly larger than the estimates from the modified 
ratio-method. Nevertheless, the estimates from DA’s using r series coefficients show a slowly decreasing trend as r 
increases. For S > 2 the highest-order DA estimates are slightly smaller than the corresponding results from the. 
modified-ratio method, but the estimate sequences show an increasing trend. If we make the reasonable assumptionEj 
that also for DA’s the dominant finite-order corrections are proportional to the amplitudes of the leading nonanalytic 
corrections to scaling, these features of the results can be simply explained by the pattern of signs and sizes of these 
amplitudes previously observed in the analysis. Taking account of these trends and performing some purely visual 
extrapolation of the DA estimate sequences, we can reconcile the DA and the analyses by the modified-ratio method 
. We shall not report in Table the DA results, but simply quote here the average of the highest-order DA data for 
a few values of S. The presence of residual trends in the sequence of DA estimates will be indicated by asymmetric 
uncertainties roughly corresponding to the range of possible extrapolations. For instance, from second-order unbiased 
DA’s, we obtain /3''“{l/2) = 0.157376|+4], /3'’“(1) = 0.224655(2), /3^='’(3/2) = 0.265640(2), /3'’='’(2) = 0.293255(2) and 

P>>^%oo) = 0.435082[+^]. 

For the same reasons as above indicated, in the sc lattice case, the values of /3/'’(<S') obtained from the DA’s are 
generally slightly larger than those suggested by the modified-ratio analysis. This upward shift is more pronounced 
in the case S = 1/2, in which we obtain /3®‘’(l/2) = 0.221665|^^p^, whereas for higher values of S the differences from 
the results by modified-ratio method are much smaller, for example: /3®“(1) = 0.312870(3), /3®'’’(3/2) = 0.368660(3) 
and /3®^(oo) = 0.601271(3). 

In conclusion, our modified-ratio method (biased with 9) and the unbiased DA estimates of the critical inverse 
temperatures on the sc and bcc lattices are consistent and compare fairly well with, but sometimes are more accurate 
than those already available in the literature for a few values of the spin-aiid^^ reported in the table. Wider 
discus^n of other estimates ^ different methods in the S' = 1/2 casaHj’ESooES, can be found in our previous 
paperEa and in recent reviewscat^J of MonteCarlo simulations and other studies of spin moi^s. It is interesting 
to mention at this point that the two most extensive simulations on the sc lattice, by a staticO and by a kineticEZi 
method, yield the estimates /3/“(l/2) = 0.2216546(10) and /3/“(l/2) = 0.2216595(15), respectively, which agree only 
within two standard deviations. 


V. MODIFIED-RATIO ESTIMATES OF THE CRITICAL EXPONENTS 


The modified-ratio methods can lead also to fairly good estimateaZa-i^the exponents 7 and v. Let us first focus 
on the calculation of the exponent 7 to recall the prescription of RefsEJ’E®!. An analogous procedure can be used for 
other exponents. For each value of S, we form the approximant sequence 


( 7 ( 5 ))„ = 1 + 


2 (s„ -f Sn- 2 ) 
i^n — 2 )^ 


where is still defined by eq.(|^) in terms of the expansion coefficients Cn{S) of xiPl S). 
Using eq. dl), we can compute the asymptotic behavior of the sequence ( 7 (S'))„ as follows 


(7('S'))„ = 7(5') 


r(7) ^(l-^^)af(5) 1 

F (7 — 6) n® n 


(26) 


(27) 
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no = 1/2, the 1/n term in eq.(p^) has a coefficient equal to (a^(S'))^ times a small positive factor. The higher-order 
corrections contain powers of 7 . As a consequence, for the Ising model, the first important correction is 0{l/n^'^^) 
and, in general, the convergence of the sequence eq.(p^) will be slower whenever the exponent under study is 3> 1 . 
This the case of 74 and, actually, we have observed that, at the presently available orders, this procedure is not 
convenient for estimating 74 , whereas a directly biased variant, to be described in the next Section, is more successful. 
On the other hand, for the calculation of the specific-heat exponent a, there are difficulties of a different nature: the 
critical singularity is very weak and the number of nonzero coefficients of the HT expansion of Cf^{/3-,S) is still too 
small. Becaus 
HT estimatesi 


we have not been able to improve by modified-ratio methods the accuracy of the current direct 
of a. 

For sufficiently large n, the sequence of approximants defined by eq.(p^ is very nearly linear on a 1/nP plot. 
Therefore, arguing like in the previous section, we are led to improve our estimates by extrapolating the odd (or even) 
subsequences linearly in 1/n®. The higher-order corrections for the exponents are expected to be more important 
than in the calculation of fdc and this reflects into a larger uncertainty of the extrapolation procedure. Just like in 
the formula for /3c, the limiting value of the approximant sequence is asymptotically approached from above, if the 
amplitude a^{S) of the leading nonanalytic confluent correction to scaling(CCS) is negative, or from below, if it is 
positive. 

In Fig. H we have plotted the approximant sequences ( 7 ®'^'^(S'))„ for several values of S between 1/2 and 00 . The 
structure of the plots is generally consistent with the pattern of signs of the CCS amplitudes already emerged from 
the study of /3/“(5'). For each sequence ( 7 ®°‘^(S'))n plotted in the figure, we have drawn as a continuous line the 
extrapolant based on the last odd pair of approximants, whereas a dashed line represents the extrapolant based on 
the previous odd pair. The small residual curvature of the approximant subsequences, which is due to the higher-order 
corrections in eq.(p7|), is made manifest in Fig. ^ by the splitting of the continuous- and the dashed-line extrapolants. 
It can also be exhibited more directly by plotting (see Fig. ^ the sequence of extrapolations of the successive odd (or 
even) pairs of approximants. 

In Table 111, for both lattices and for several values of the spin, we have reported the numerical values of the 
extrapolated exponents with an error corresponding to a small multiple of the difference between the continuous and 
the dashed extrapolations. We have also reported the derivatives of these estimates with respect to 6, computed at 
the reference value = 0.504. For comparison, the same Table also shows the exponent estimates obtained from 
DA’s, while the results obtained in other recent numerical studies using MonteCarlo methods, by shorter HT series. 


or in the RG approach, will be further discussed in Sec. VHI and IX and are collected in the next TabletV, 

From Figs. ^ and it is clear that not only a®°^(5'), but also the amplitudes of the main subleading CCS change 
sign as S varies between 1 and 2. This very favorable circumstance, which can also be confirmed numerically, for 
example by fitting the approximant subsequences to the simple asymptotic form 7 -|-Ci(S')/n® + C 2 {S)/n^~^^, makes us 
very confident about the accuracy of the exponent estimates presented below. For each value of 5, a simple monotonic 
behavior appears to have set in, since as shown in Fig. the subleading asymptotic correction in eq.(p^ generally 
works in the expected “right” direction. Namely, it tends to lower the extrapolated exponent values obtained from 
the decreasing approximant subsequences for spin S' = 1/2 and S = 1 , while it tends to raise the extrapolated values 
obtained from the increasing subsequences for S >2. Only in the S = 3/2 case, in which both the amplitudes of the 
leading and of the subleading correction have the smallest absolute value, the approximant sequence is very slowly 
increasing and the sequence of extrapolated exponents is very slowly decreasing. Thus we can expect that, as the 
number of available coefficients grows large, the range of variation with respect to S of the extrapolated estimates of 
7 ®““(S') will continue to shrink, further improving the verification of the universality of the exponent with regard to the 
spin. More precisely, assuming that the general features of the behavior we have described persist as the order of the 
series increases, the successive extrapolations of the sequences ( 7 ®°°(l/ 2 ))„ and ( 7 ®°“(l))n should provide decreasing 
sequences of upper bounds, while those of the sequences ( 7 ®“( 2 ))„, ( 7 ®“( 5 / 2 ))„ etc. should give increasing sequences 
of lower bounds for 7 . 

At the present order of expansion, the exponent estimates obtained by our extrapolation prescription, range orderly 
from 1.23742, for S = 1/2, to 1.23684 for S = 00 . Therefore, if we now assume that universality is valid, in particular 
that 7 is independent of S, the previous remarks suggest to take simply the average 7 = 1.2371(4) of these extrema as a 
first rough approximation of the exponent with an uncertainty corresponding to the half-width of the range of variation. 
We can further refine this estimate observing that, for values of the spin between 1 and 2, both the leading and the 
main subleading CCS are very small, as it appears observing that the exponent approximant sequences have very small 
slopes, clearly positive for 5* = 1 and negative for S = 3/2 and S = 2. Moreover the extrapolated exponent estimates 
are fairly insensitive to the bias value of 9 (for instance, we have d^^^^{l)/d9 ~ 0.003, 57®“(3/2)/90 « —0.0005 and 
dj^^^{2)/d9 « —0.002 at 0 = 0”®-^). Since ( 7 ®''^''^(l))rt and ( 7 ®“( 2 ))„ are very close, a better estimate for 7 should 
lie in between. The extrapolation of the last odd pair of terms in the sequence ( 7 ®''^''^(l))r!, yields 1.23730, whereas 
for the sequence ( 7 ®“( 2 ))„ it leads to 1.23699, and therefore the rough estimate given above can be improved to 
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7 = 1.23715(15). Consideration also of the sequence ( 7 ^'^‘^( 3 / 2 ))„ suggests that we take 7 = 1.2371(1) as our final 
best estimate. . 

A closely related procedure was proposed long ago in RefsBa. These authors analyzed but never published, extensive 
two-variable series in power of (3 and of a continuous Ising spin variable (made available by B.G. Nickel), using partial- 
differential approximants methods which indicated an “effective fixed point” around S = 3/2. Within the precision 
of the present calculations, the simple prescription of taking the average of the extrapolations of ( 7 ^““(l))n and of 
( 7 ^“'=( 2 ))„, or the extrapolation of ( 7 ^“( 3 / 2 ))„ as the best approximation of 7 , should be equally (or perhaps more) 
effective since it also takes advantage of the vanishing of the main subleading correction. 

Since it is not difficult to show that the leading CCS for any observable must also vanish for the same value of A, the 
same prescription can be used for extracting the best value of v from the approximant sequences formed by 

the series coefficients of ^^(/3; S) and shown in Fig. The following Fig. shows the sequence of the extrapolations 
of the successive odd pairs of approximants. The slower convergence of the approximants to the correlation-length 
exponent should not be surprising, simply because /i 2 (/ 3 ; 5 ') enters into the definition of At the present orders of 
expansion, the behavior of the sequence of the extrapolated exponent values is clearly not yet asymptotic for 5" = 1/2 
and S' = 1, while it is much smoother and shows a slowly increasing trend for 3/2 < S < 3 and a slowly decreasing 
trend for S > 3. Thus arguing as before, we can conclude that v = 0.6299(2). 

If we bias the extrapolation procedure with a larger value of the range of variation of with S will be 

expanded, to an extent that can be easily figured out from the data reported in Table 111, but the estimated central 
value of 7 ^““ will be practically unchanged. For instance, if we adopt the significantly larger value 0”®-^ = 0.54, we find 


7 ''“( 1 / 2 ) = 1.23782 and 7 ''“(oo) = 1.23661, whereas 7 ''“( 1 ) = 1.23741, 7 ''“=( 3 / 2 ) = 1.23708 and 7 ''“=( 2 ) = 1.23691 
are changed to a smaller extent. Averaging and yields 7 = 1.23716(25) and consideration also 

of ( 7 ^''^''^( 3 / 2 ))„ leads to essentially the same final estimate as the one obtained for 0”®-^ = 0.504. In the case of 
the exponent the estimated central value is slightly lowered to 0.6298, well within the error bars of our previous 
estimate. — 

Using the Fisher scaling relationOBil, the exponent rj describing the large distance falloff of the two-spin correlation- 
function at the critical temperature can be estimated 77 ^“ = 2 — = 0.0360(8). 

In Fig. 1^ and || we have shown the results of the analogous procedure of extrapolation for ( 7 ®°(S'))„ and {u^^{S))n- 
The main features are similar to the bcc case, except, unfortunately, for the sign pattern of the amplitudes of the leading 
CCS, all of which now appear to be negative, consistently with the study of /3/“(<S') by the modified-ratio method. In 
complete analogy with the bcc case, for both exponents 7 and ly, the residual curvature of the approximant sequences 
tends to correct the extrapolations in the direction expected in order that universality be realized. The accuracy of 
the exponent estimates is notably smaller, due to the anticipated slower convergence of the sc lattice series coefiicients 
to their asymptotic form and perhaps to the presence of larger subleading corrections in eq.(27). The values of the 
exponents 7 and iz obtained from the sc series data are consistent with those from the bcc data, but they are affected 
by significantly larger uncertainties: we can roughly estimate 7 = 1.2368(10) and v = 0.6285(20). 

The numerical process achieved in this study is best apraeciated by comparing our Figs. ^ and|| with the analogous 
Figs. 1 and 2 of Ref.Efl. We should first observe that in Ref.cJ a straightforward extrapolation linear in 1/n was implied 
for the sequences and (i/^““(5'))„. Due to this choice of the plotting variable and to the smaller extension of 

the bcc series available two decades ago, the “relative maximal spreads” with respect to the spin S, of the extrapolated 
exponent values are !^| 2 / 2 )+ 7 (m) ~ ^ ~ ^ 10“^, respectively. In our study of the same 

lattice, the corresponding figures are smaller by nearly one order of magnitude, namely the relative spread is now 
« 2.3 X lO”"* in the case of 7 and « 1.4 x 10“^ for zz. The values of these spreads can be taken as rough accuracy 
limits for the verification of the universality with respect to 5, which is thereby convincingly corroborated by the new 
analysis. 

We close this discussion with a few remarks. Our extension of the series to order has been crucial in showing 
that, in the bcc lattice case, the asymptotic structure of the HT expansion coefficients is already well stabilized, 
since the last six or seven modified-ratio method approximants of the critical inverse temperature or of the critical 
exponents show remarkably regular trends. Also in the sc lattice case, there are indications from the last three or 
four approximants obtained by the same method, that a similar trend is setting in, but clearly the convergence is not 
as fast as for the bcc lattice. 

Some numerical experimentation with model series suggests, also for the exponent analysis, that our error estimates 
are reasonable and quite conservative. 

For both lattices, the CCS amplitudes can be estimated from the slopes of the exponent approximant sequences, 
as will be further discussed in Sec. IX. 

In conclusion, this simple modified-ratio approach confirms accurately the universality of 7 and v with respect to 
the magnitude of S and to the lattice structure and, conversely, assuming universality and using the bcc series data, 
it yields very accurate estimates for these exponents. 
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VI. BIASED MODIFIED-RATIO METHOD FOR THE EXPONENTS 


In Ref.Ea, J. Zinn-Justin proposed also a more direct modified-ratio procedure for biasing the exponent estimates 
with the value of 6, in order to eliminate or strongly reduce the influence of the leading confluent corrections to scaling. 
The prescription involves the quantities 


Sn — {Sn + Sn-l)/2 

(28) 

and 


K = (1(.U - 

(29) 

In terms of bn the following approximant sequence can be formed 


1 1 {bn ~ bn-2) 

(30) 

If we make the simplifying assumption that 0 is exactly 1/2, also the correction terms 0(l/n^® 
by this prescription, along with the regular correction 0{l/n) and therefore 

) will be eliminated 

(7(*S'))n -7(>S') + 0(^3/2) 

(31) 


By the remarks made at the beginning of the preceding section, these are not decisive improvements in the calculation 
of 7 and and indeed, both for the sc and the bcc lattice, the results obtained by this procedure are consistent with 
but not more accurate than those of our previous analysis by modified-ratio methods. See for example Figs. || and 
p^ , where, for convenience, we have plotted {^{S))n vs. l/n^+® rather than vs. because the plots of the 

approximants appear to be more nearly linear (although with somewhat large corrections) with respect to former 
than to the latter variable. 

On the other hand, this biased variant of the modified-ratio method, is more successful in the analysis of the 
expansions that we have computed for 74 . The sequences of biased approximants {lf{S))n for the exponent of 
x'f (/3, S)^ are shown in Fig. |l^ for the bcc lattice and in Fig. p^for the sc lattice. In order to avoid confusing the plots, 
we have indicated only the extrapolations, linear in based on the last odd pair of terms {( 74 ^(S')) 2 i, {xf (<S')) 23 } 

in the approximant sequences. When the spin S varies between 1/2 and 00 , the extrapolated values of 7 |'^( 5 ') range 
from 4.366 for S = 1/2, to 4.372 for S = 00 . Similarly, in the case of the bcc lattice the values of 74 “(S') vary between 
4.369 and 4.375. We have reported in Table ^ the results obtained by this method for several values of S. From the 
sc lattice data we can conclude that 7!'^ = 7-I- 2A = 4.369(8) and from the bcc lattice data 7^'^'^ = 7-I- 2A = 4.372(8). 

The accuracy in the verification of the validity of hyperscaling is often characterized quantitatively by quoting the 
value of the right-hand side of eq.(|l 8 |): from our estimates we have: 7 J- 3z^ — 2A = 27 J- 3;/ — 74 = —0.0099(160) in 
the sc lattice case, and analogously j + Sv — 2A = —0.0081(88) for the bcc lattice. 

These results give strong support to the validity of the hyperscaling relation and of the universality of 74 with 
respect both to the lattice structure and to the value of S. 


VH. RATIO ESTIMATES FOR THE EXPONENT OF THE LEADING CONFLUENT SINGULARITY 

Assuming that 0 is universal, the simplest prescription for estimating this exponent is based on the series with 
coefficients 


qn{Si,S2) 


Cn{SMS2) 

CniS2)dn{Si) 


(32) 


for n > 0 and, of course, Si ^ S 2 - Here c„(S') are the coefficients of the susceptibility and d„(S') the coefficients of 
the second correlation moment for spin S. From eq. (H) we can observe that, for large n 


qniSi,S 2 ) = A(5i, S 2 )[l + + 0{l/n)] 


(33) 


therefore 
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qn{Sl,S2) 
qn+2iSl, S2) 


1 + 


A' 


+ 0 (l/n^) 


(34) 


so that 


(0(Sl,S2))n 



rn(Sl,S2)-l 
‘fn+2(S'l, S 2 ) — 1 


l )-2 


(35) 


is an approximant sequence for 0. In the bcc lattice case, if we choose Si = L(2, S 2 > 2, and extrapolate only the even 
(or, equivalently, the odd) subsequences linearly in l/n^+®, we obtain Fig. |^. The results indicate very suggestively 
that 9 — 0.50ti, independently of the value of 82 - 


VIII. AN ANALYSIS OF THE EXPONENTS BY DIFFERENTIAL APPROXIMANTS 


The modified-ratio methods employed in the last sections have proved successful and suggestive both for the 
determination of the critical temperatures and for the calculation of the exponents 7 , v and 74 . Let us now turn 
to the more traditional differential approximants (DA) based procedures of series analysis after recalling that their 
main difficulty is the necessity of some further extrapolation with respect to the order of the series used, which is not 
straightforward, due mainly to the lack of simple estimates for the finite-order corrections and to the spread of the 
various DA estimates at a given order of approximation. This fact also hampers the assessment of the errors, which 
can be realistic if not only they reflect the spread of the values of the highest-order approximants, but also allow for the 
possible residual trends. In this respect, the modified-ratio methods might be easier to use, as we have suggested in the 
previous sections. We have already discussed in Sect. IV the DA estimates for the critical points. For measuring the 
exponents, we have preferred series analyses using the inhomogeneous first- and second-order DA’s biased with /3c (or 
in some cases with ±/3c), or sometimes the simplified inhomogeneous first-order differential approximants defined ino, 
in which we have fixed also the correction-to-scaling exponent 9. The extrapolations of the results from the biased 
DA’s and from the simplified DA’s may be performed with a smaller uncertainty, because the spread of the estimates 
tends to be narrower than for unbiased approximants. Moreover, in order to understand, at least qualitatively, how 
the estimates on a given lattice depend on the spin and to improve them, it will be sufficient to assume that the leading 
finite-order corrections-arje.^roportional to the amplitudes of the leading nonanalytic corrections to scaling(CCS). 

A simpler approacht23’t23 similar to the simplified DA’s consists in forming the conventional Fade approximants 
after subjecting the series to the biased variable change w^(S') = 1 — t^( 5')® in order to regularize the leading 
CCS. The results obtained either by simplified DA’s or by Fade’ approximants in the variable w'^{S) are sometimes 
numerically comparable, but the latter are generally affected by a larger uncertainty. 

We have computed also the ffi^tivej-eaponents, introduced long ago in Ref.tHd and more recently reconsidered and 


systematically studied in Refs. I 


L0f,109 


for the susceptibility 




d\ogx*ifi;S) 

dlogT#(^S) 


= j{S) - 9a*{S)T*{Sf + OiT*{S)) 


(36) 


for the correlation length 

4//3; S) ^ -(1 - = 12 ( 8 ) - (5)r#(5)^ + 0(r#(5)) 

and for the second field derivative of the susceptibility 


itfM.s) 


dlogxf (/3; 8) 
dlogrib i^S) 


74 (^) - 9at{S)T*{Sf + 0{t*{S)). 


(37) 


(38) 


The critical exponents 7 , v and 74 are estimated by extrapolating the effective exponents to the critical singularity. 
Of course, the factor (1 — t(S')) in eq.(^ is introduced only to compensate for the singularity of at /3 = 0 

and is unimportant at the critical point. 

It is interesting to plot the effective exponents over a wide vicinity of /3^(5'), not only to gain informations on 
the leading correction amplitudes af'(S'), a,f{S) and af{S) through eqs.(|^)-(^, by examining whether and how 
fast they approach the critical limit from above or from below, but also simply in order to display the varietyE. 
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of preasymptotic critical behaviors which can occur within the same universality class. The parametrizations of the 
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approachestiilliii 


to RG, must confront 


approach to the critical behavior, proposed within various field-theoretical 
also with this phenomenology. 

In Fig. lU and ^ we have shown the highest-order simplified DA’s of the effective exponents 7 gj^(S') and, respec¬ 
tively, ly^ff (S) for spin S = l/2,l,...oo, over wide ranges of inverse temperatures. In order to make the curves 
conveniently comparable for all values of the spin, we have plotted the effective exponents against the variable 
(.j.bcc(^))e _ rpjjg leading CCS is revealed by the slope of the plots near = 0. 

While the simplified DA’s are quite sufficient to give a general view of the behavior of the effective exponents, more 
accurate results for the exponent 7 are obtained extrapolating the effective exponent expansions by second-order 
inhomogeneous DA’s biased with f3c. The estimates thus obtained for 7 ( 5 ') range from 7 ^°°(l/ 2 ) = 1.2385(6) to 


7 '>“(oo) = 1.2367(5). 


averaging 7 ^°°( 1 ) and 7 ^“( 2 ) and taking into account also the value of 7 ^°°( 3 / 2 ). It agrees well with the estimate 
by modified-ratio methods. The corresponding results for the correlation-length exponent range from z/^°°(l/2) = 
0.6314(20) to = 0.6294(5) and our best estimate is v = 0.6301(4). 

In the sc lattice case, the analogous (but less well converged) plots for and for 5), obtained by 

simplified DA’s, are shown in Figs. and |^. This analysis also confirms that, on the sc lattice, the amplitudes of 
the leading CCS do not have a dependence on S similar to the bcc case, but remain negative for all values of the spin. 
The estimates of 7 '*°(S') and of obtained by second-order biased DA’s are also reported in Tableffl . 

The simplified-DA analysis of the effective exponent 74 );^^(S') of X 4 “(/ 3 ;S) yields estimates of ranging 

between 4.3647 and 4.3653. It also indicates that a 4 “(l/ 2 ) and 04 ^( 1 ) are negative, whereas 04 ^( 5 ) is positive for 
S > 3/2. On the sc lattice, the corresponding estimates of 74 °(S) vary between 4.363 and 4.373 and a 4 °(S) > 0 for 
all S. We can conclude that 74 = 4.366(2), independently of the spin and the lattice, and in good agreement with 
the results of the analysis by modified-ratio methods. The accuracy in the verification of hyperscaling is now slightly 
improved with respect to the biased modified-ratio methods of sect. VI, since we have 7 -I- 3z^ — 2A = —0.0021(28). 

As shown in Figs. and |^, the pattern of signs for the confluent amplitudes of consistent with the 

corresponding results for x#(/3;5'), as it must, since the ratios af {S)/a^{S) are expected to be universal. 


They are reported in Table 

,bcc 


III 


Our best DA estimate 7 = 1.2373(4) is obtained simply by 

..beef 


The exponent xt (S) can also be evaluated extrapolating the effective exponents by inhomogeneous second-order 
DA’s biased with On the bcc lattice, our results, which appear in Table 0, range between 4.376(8) for ^ = 1/2 

and 4.3631(4) for S = 00 . In particular we find 74 ^( 1 ) = 4.3666(10), 74 ''’''’( 3 / 2 ) = 4.3638(10) and 74 ^( 2 ) = 4.3629(10). 
By the same arguments used in the modified-ratio method analysis of the bcc lattice series, the best value for the 
exponent 74 should lie between the estimates for S' = 1 and S = 2. This leads to our final estimate 74 = 4.3647(20), 
which is even more accurately consistent with hyperscaling, since for the bcc lattice data we have 27 -|- 3:/ — 74 = 
-0.0008(28). 

Finally, it is worth to mention briefly that the ratio iz/y can be determined, to a good precision, also studying the 
log-derivative ratios Dlog{xi)/Dlog{x) and Dlog{g, 2 /[})/Dlog{x), either by DA’s biased in Pc or by simplified DA’s 
biased in Pc and 9. The values thus obtained from the bcc lattice expansions (except for S = 1/2) fall within the error 
bars of our best result iz/y = 0.5092(2) from modified-ratio methods. The accuracy of the estimates can be further 
improved by focusing on the bcc lattice case and arguing as usual that the best value of 1^/7 is simply an average of 
the estimates for S' = 1 and S = 2. We thus arrive to the value vjx = 0.5091(1). The estimates of this ratio obtained 
from the sc series lie within twice the expected error bars for S > 2, but are slightly worse for smaller values of S. In 
the bcc lattice case, also the DA estimate of 7 obtained from the analysis of the ratio of the log-derivatives of x! and 
X, whose value at the critical point is 1 -|- I/ 7 , agrees very closely with our best results by modified-ratio methods. 

We have also examined the term-by-term divided series 


Q{x-S) 


E 


Pr{S) 


(39) 


where er{S) are the expansion coefficients of X 4 (/ 3 ;S') and Pr{S) those of x^{P\S). Using eq. ( 11 , it is easily shown 
that the auxiliary function Q{x', S) has a critical point at a; = 1 with exponent Zv 1, if hyperscaling holds. A 
second-order biased DA analysis of the effective exponent yields the estimate v = 0.6300(4) independently of S and 
in complete agreement with hyperscaling. 

In conclusion, we have observed that if the sequences of modified-ratio method approximants are carefully extrapo¬ 
lated using as bias the value of 9 derived by RG methods, the estimates of the exponents obtained by the modified-ratio 
method, for all S and on both lattices under consideration, show a good agreement with the results from DA’s biased 
only with Pc or simplified DA’s biased with both Pc and 9. The efose consistency between the critical parameter 
estimates obtained by a number of different procedures adds further confidence that the HT series have now reached 
a fairly safe extension and that we are not being misled by accidental apparent convergence, so that the uncertainties 
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of the HT estimates can be significantly reduced. The small residual dependence of the exponent estimates on the 
spin S and on the lattice structure can be confidently used to characterize how accurately universality is respected. 


IX. A COMPARISON WITH OTHER EXPONENT ESTIMATES 


The agreement of our HT estimates with the values 7 = 1.2396(13) and ly = 0.6304(13), obtained 1 
RG by Borel-summation of coupling constant seventh-order fixed-dimension perturbative expansion 


i^ontext of 
_£;.r-pr with 

the values 7 = 1.2380(50) and iz = 0.6305(25), obtained by Borel-summation of the fifth-order e— expansionES^, is still 
reasonable. The values 7 = 1.2378(6) -I- 0.18((jfr- — 1-40) and ly = 0.6301(5) -I- Q.12{gr — 1.40) proposed in Ref.i£ on the 
basis of a slightly different resummation of the fixed-dimension perturbative RG expansion, are perhaps even closer. 
At the presently available orders of HT expansion, our series estimates prefer central values for 7 and v which are only 
slightly lower. It is appropriate to meptionrtiiat very similar central estimates, though with largef|riiiicertainties, were 
already obtained quite some time agcc§ETE3 from bcc latoce series of order by the methodc^Ea of Ghen, Fisher, 
Nickel-apd Zinn-Justin. For instance, the analysis of Ref.cJ yielded 7 = 1.237(2) andr^p^^300(15). A more recent 
studyllia of HT series through 0(/3^°) for the sc lattice, along the same lines as in Refs.E3Ea’E£l, indicates 7 = 1.2371(4) 
and V = 0.63002p|3^ All these results are also quoted for comparison in our Table 0 

The techniquecjcj of focusing the analysis on some particular model in the Ising universality class with negligible 
amplitudes of the lea^g confluent corrections to scaling(CCS) was advantageously adapted also to MonteGarlo 
simi^tions in Refs.Ou, which report 7 = 1.2372(17) and v = 0.6303(6). This procedure was further improved in 
Ref.H, in which it led to the estimates v = 0.6296(7) and 77 = 0.0358(9), implying 7 = 1.2367(20). Even lower central 
estimates of the exponents, namely 7 = 1.2353(25) and v = 0.6294(10) have been measured in a npse conventional 
MonteGarlo simulation of the spin-1/2 Ising model supplemented by a finite-size scaling analysisiil which allows 
also for the corrections to scaling. It is tempting to conjecture that, for S' = 1/2 on the sc lattice, our results from 
the extrapolations of the modified-ratio approximants and the best finite-size-scaling analyses of the MonteGarlo 
simulations on the largest accessible lattices are subject to errors of the same nature. This would explain the rather 
small central values of the quoted MonteGarlo estimates of 7 and indicate the need of simulations of a significantly 
larger scale in order to obtain from spin- 1/2 systems on the sc lattice exponent values in closer agreement with our 
bcc-series estimates. 

Let us now comment briefly on the existing results-|fpC| 74. We recall that the validity of the hyperscaling relation 
eq. (H) for the spin-1/2 Ising model was questionedETTa on the basis of an analysis of 10-12 term series on the sc, 
the bcc and the fee lattices, yielding the estimate 7 -I- 3 z 7 — 2A = 0.038(12). This result was at the time interpreted 
as an indication of a small, but clear, failure of hyperscaling. As already mentioned in the introduction, the problem 
was convincingly settled only the HT series for y and on the bcc lattice, extendedE3 up to order 21, were 

analyzed with careful allowanceEJTa for the GGS and indica|te(^^ insufficient accuracy of the “classical” HT estimates 
7 = 1.250(3), V = 0.638(2) and a-p 1/8 generally acceptecEllSEa until then. 

For general spin, a single studjcd of X 4 , performed with series 0(/3^^) on the fee lattice, can be found in the literature. 
The log-derivative of the series Q{x] S'), defined by eq. (^), was examined by PA’s. The analysis produced estimates 
of the exponent 2A — 7 ranging from 1.887(1), for S = 1/2, to 1.893(1), for S = 9/2. Of course, if hyperscaling is 
valid 2A — 7 = 3 z 7 . Thus the^fu^l estimate 2A — 7 = 1.890(3) indicated the absence of hyperscaling violations of 
the size predicted from RefsoLj, provided that the central value v = 0.630 suggested by the RG, rather than the 

“classical” HT estimate v = 0.638(2), was adopted. _ 

The expansion of X 4 on thOpSC lattice, for S = 1/2, was at the time already availableE/j up to order 17, but it 
was analyzed only later in Ref.c3. It yielded the estimate 74 = 4.370(5) still confirming the validity of hyperscaling, 
provided that the revised values obtained from the RG in those years, weije assumed for 7 and v. The series on the 
bcc lattice, for S = 1/2, was extended to the same order only much laterE^ and its analysis also cpefiipied the above 
conclusion. Further support of these results came also from various more recent MonteGarlo testsdJilia. 

It should be stressed that the finite-order effects are stronger in the calculation of X 4 (and of the related quantities 
like gr)— than in the calculation of quantities defined in terms of two-spin correlations, as we have already remarked 
in Ref.u, and therefore that the accuracy of the results is correspondingly smaller. Our comparison with previous 
studies, shows, however, that we have achieved some improvement not only in the precision of the estimates of 7 and 
V, but mainly of 74 , by taking advantage of our significantly extended expansions of X 4 - 


14 








X. ESTIMATES OF CRITICAL AMPLITUDES 


For proper reference and for comparison with the earlier studies, we have reported in Table VI a set of updated 
estimates for the critical amplitudes of x^(/3; S), of S), of S) and of xt (/^! defined by eqs.(||)-(|2|). 

We have evaluated the critical amplitude (7^(5') of x^{/3;S') as follows. We have adopted as a bias the value 
7 = 1.2371 and our estimates of l3f{S) to compute the HT series of the “effective amplitude” 


C*{P; S) = T^x*{P; S) C*{S) (l + a*{S)T*{Sf + ... + b*{S)T*{S) 


m 


The amplitude C'^(5') is then estimated by extrapolating the effective amplitude S) to the critical point. The 

extrapolation has been performed either by first-order inhomogeneous simplified DA’s biased with fjf (S) and with 6 in 
order to allow for the confluent corrections to scaling or, more traditionally, by using second-order inhomogeneous DA’s 
biased with ±Pf{S). Since these two procedures yield fully consistent estimates, we have reported in Table VI only 
the results obtained by the usual differential approximants, which do not need to be biased also with 9. By the same 
procedure, we have also studied the effective amplitudes for the correlation length and for Xi in order to evaluate the 
corresponding critical amplitudes (/'^(S'))^ = limi- ^o+T^‘'^*{P;S)^ and CfiS) = - lim,^o+(/3; 5). The 
results are reported in the same Table. 

The above mentioned difficulties in the analysis of the critical behavior of the specific heat, also result into larger 
errors of the critical amplitude A^{S). Therefore it seems more convenient to compute this quantity from the second 
derivative of S), which presents a sharper singularity. 

Other estimates for some of the mentioned critical amplitudes, obtained from shorter HT series apJ. rmder-shghtly 
different biasing assumptions, or by other numerical methods, can also be found in earlier studiesE3’Bli£3ii£3. For 
instance, from Ref.ll23 we have cited in Table VI the estimates of A®‘^(l/2), C'®'^(l/2) and /®‘^(l/2), obtained from 
series 0(/3^’^), and respectively, under the assumptions a = 0.104, 7 = 1.237, ly = 0.6325 and 

/3®°(l/2) = 0.221620. From the same reference, we have also reported the estimates of A*'“(l/2) obtained from a 
series 0(/3^'^), and of C'*'“(l/2) and /*'“(l/2), obtained from series 0{/3^^), by assuming /3j“(l/2) = 0.157362 and 
the same values as above for a, 7 and ly. _ 

Under various assumptions on the value of a, estimates of A®°(l/2) were derived in Ref.EHj from a simulation in 
which the energy and the specific heat were measured. By straightforward interpolation, we can conclude that, for 
a = 0.11, these, data would imply the estimate A'*°(l/2) = 1.368(7), in reasonable agreement with ours. 

From Ref.tHJ, we have quoted estimates of C'|°(l/2) and of C4‘^^(l/2) obtained from series 0(j3^^) and 0(/3^^), 
respectively, assuming 74 = 4.375 and /3f (1/2) = 0.221630Lm), /3/“(l/2) = 0.157368(7). 

In the same Table, we have also reported the estimates^ from series for C'^“(l/2), and C'^°°(2) 

assuming 7 = 1.237 and the estimates of /^'’’^(l/2), /^^‘^(l) and /^“(2), obtained assuming ly = 0.6297, ly = 0.6298, 
1 / = 0.6300jespectively, together with the values of quoted in the same Reference and reported in Table ||. 

In RefJlH3 the values of C'®‘^(l/2) and of f®‘’(l/2) have been computed by a MonteCarlo method assuming 7 = 1.237, 
1 / = 0.628 and /3f (1/2) = 0.22165. 

This brief review of some existing results shows how sensitively the estimates of the critical amplitudes depend on 
the bias values chosen for the critical exponents and temperatures, and, of course, on the length of the series. If we 
also allow properly for these sources of uncertainty, which generally are not included in the error bars quoted in the 
literature, many of the cited estimates can be considered essentially compatible among themselves and with ours. 


XI. ESTIMATES OF THE RENORMALIZED COUPLING 

The value of the hyper-universal renormalized coupling constant gr can be obtained from our estimates of the critical 
amplitudes. Alternatively, without biasing the computation also with the critical exponents, gr can be computed, 
with a smaller uncertainty by extrapolating to critical point the expansion of the auxiliary function 

y*ip-,S) = {g*ip-,S))-^/^ (41) 

or of the function 

S) = {(3/(3*{S)f/^g*{p-, S). (42) 

Unlike the effective coupling gf{(3]S), both y^(/3;5') and z'^{f3]S) are regular analytic at /3 = 0, so that they can 
be expanded in powers of (3 and extrapolated to the critical point by Fade approximants, DA’s or simplified DA’s. 
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Due to the finite extension of the series, the numerical estimates of gr derived from eq. (© or from eq.(^^ are of 
course (very) slightly different. In order to allow for the expected leading confluent corrections to scaling(CCS), in 
our calculations we have used first- and second-order DA’s biased in Pf{S) or simplified DA’s biased with Pf{S) and 
with the confluent exponent 9. In Fig. |2^ we have plotted vs. the effective coupling S) as obtained 

from the function S) for various values of the spin S. The curves, computed in the simplest way by simplified 

DA’s, show the strong influence of the CCS nearby the critical point and indicate that is independent of S 

within a very good approximation. Comparison with Fig. which shows the effective coupling S) plotted vs. 

similarly indicates that the renormalized coupling 5 ^° « « 1.41, is independent not only of S, but also 

of the lattice structure. Using eq. dil), we can infer from Fig. ^ that the amplitudes of the CCS are generally large 

Fig we 


IHJ in the 


and, more precisely, that ag“(S') > 0 for S' = 1/2 and 1, whereas ag“(S) < 0 for S > 2. Analogously, froipJ 
can conclude that ag‘^{S) > 0 for all S. These qualitative conclusions are consistent with RG estimates^ 
fixed-dimension approach indicating that agjay- lies in a range from —3. to « —2.. 

In order to reach higher precision in the calculation of gr, we have preferred to use first- and second-order DA’s 
biased with /3c. In the bcc lattice case, we notice that, for S = 1/2 and 1, approximants which use an increasing 
number of coefficients show a residual slowly decreasing trend, while, for S > 2, they show an increasing trend. We 
shall indicate by asymmetric uncertainties these features of the approximant sequences. Again arguing like for the 
critical exponents 7 , v, and 74 , we can expect that the most reliable estimate of gr obtained from the bcc lattice 
series will be nearly equal to the average of g/“(l) and (;/°°(2) or to the value g/“(3/2). Thus w 
gr = 1.404(3). The central value of our updated estimate is slightly lower than our previous resultl^ 
based on a series 0(/3^^) for the bcc lattice in the S' = 1/2 case, in which the convergence is slo; 


est. However, our 

revised result is slightly closer to the value gr = 1.400 advocated by D.B.Murray a^id^B.G.Nickels in the context of 
the RG fixed-dip 
method of Refs.c 


conclude that 
gr = 1.407(6), 


ion approach, and is compatible with the more recent HT resultHi^ gr = 1.402(2) obtained by the 
from sc lattice series exten ding to order /3^®. Our numerical estimates of gf{S), for several values 
of S and on both lattices, are reported in Table VII. Notice that in the bcc case the DA estimates are larger or smaller 
(in the sc case generally larger) than the expected best value, consistently with the sig ns of th e CCS amplitudes. In 
the same Table we have quoted for comparison also some recent MonteCarlo estimatesIliSlliS of gr 
results from HT and RG methods. 


as well as other 


XII. ESTIMATES OF A? 


The combination of critical amplitudes , defined by eq. (|^), can be computed either from the estimates of the 
critical amplitudes and of the exponent a = 2 — ?>v, or, more directly, by extrapolating to the critical point 

the expansion of the auxiliary function 


F{(3;S) = 


-qv^f3r{S)d'^G{(3-S) ( 


2 vl 3 


d/32 [02 


/ 33/2 d{l/aP;S)) ^ 

(S) d(3 ) 


-3 


(43) 


since the validity of the hyperscaling relation eq.(20) implies F{Pc{S) — 0; S) = (i?^)^ + 0(l/n®). We have assumed 
= 0.6299 and evaluated F(/3; S) by first-order d iffere ntial approximants biased with The estimates of 


obtained by this prescription are shown in Table VHI. Within a fair approximation, they are independent of S and 
of the lattice structure and compatible with the estimates obtained combining the amplitudes reported in Table|V4| . 

Our final estimate is R^ = 0.2668(5). This result is slightly smaller than that the value R^ = 0.272(4) reported in 
our previo us stu d>@ of the single S' = 1/2 case employing shorter series. 

we have also shown values of Rt obta ined via RG, either by the fixed-dimension perturbative expansion 


In Table VIH 


to fifth ordeillHa or by the e-expansion to second ordeiUHS. We have also quoted the estimate R^ = 0.270(4) that would 
be obtained from the MonteGarlo measures of Ref.tla assuming a = 0.11 and the recent HT result R^ = 0.2659(4) 
taken from Ref.0 as a representative of various nearly equal central estimates from studiesiliiS performed at 
different times, with different techniques, under different assumptions on the values of v and a and using series of 
different extensions. The discrepancy from our estimate should probably be taken as an indication of the remaining 
difficulty of accurately evaluating the specific heat amplitude. 

Finally, it is worth while to quote t wo r ecent very accurate measurements on binary mixtures: 37^ = 0.2841181 
performed in a microgravity experimentll2Zl, and Rt = 0.265(6) obtained in a conventional environment 
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XIII. ESTIMATING THE RATIOS OF CONFLUENT-SINGULARITY AMPLITUDES 


From the extended series presented here, we have also tried to evaluate the universal ratio af(S)/a^(S). We recall 
that, for the bcc lattice, our analysis of x and ^ by modified-ratio methods had shown that, as the spin S varies between 
1 and 2 , the leading correction amplitudes a^{S) and a^{S) vary from small negative values to small positive values, 
whereas in the sc lattice case no change of the sign of the confluent amplitudes is observed. As we have emphasized, 
some knowledge of these amplitudes is necessary to understand how the various numerical estimates obtained for each 
value of S approach the true values of the universal critical parameters. A simple prescription to compute accurately 
the universal quantities consists in using series on the bcc lattice with spins between 1 and 2 , for which the amplitudes 
of the leading confluent corrections to scaling(CCS) are very small. Conversely, the numerical methods to evaluate 
the amplitudes and the exponent of the leading CCS, can be expected to work with fair accuracy only when the 
confluent corrections are not too small. For S = 1/2 the size of the leading CCS is largest, but unfortunately also the 
higher-order corrections seem to be still important, as shown by the steep behavior of the extrapolated sequences in 
Fig. 1^. Therefore, the most reliable results are likely to come from the bcc series for S > 2. 

We have obtained reasonably accurate estimates _of the CCS amplitudes for x simply by fitting ^ asymptotic 

under the 
assuming 


form 7 -I- ci{S)/n^ -I- C 2 (S')/n^+®, suggested by eq. (izD> to the exponent approximant sequences eq.( _ 

assumptions that 7 = 1.2371 and 0 = 0.504. A similar procedure can be repeated in the case of 


V = 0.6299. The estimates of the CCS amplitudes and of the universal ratios a|“(«S')/a^°'^(S') thus obtained, are 


shown in Table tX. The values of the amplitudes obtained by this procedure are generally consistent within 

a few percent, and those of within 10 — 20% in the worst cases, with those evaluated by simplified DA’s of 

the log-derivatives of x and biased with /3c and 9. Moreover, for S' > 2, the estimates of a^“(S) obtained from the 
modified-ratio approximants for 7 are consistent within a few percent with those obtained from the corresponding 
approximants for /3c. As expected, for smaller values of S, we have consistency only within 20 — 40% because the rate 
of convergence of the series is slower and/or the subleading confluent corrections are more important. 

The ratios aj°‘^(S)/a^°^(S) appear to be approximately independent of the spin S, as they should, and suggest the 
final estimate a^/a-^ = 0.76(6). The error includes also the uncertainties of the bias parameters 7 , ly and 9. In the sc 
lattice case, the same analysis leads to amplitude ratios which show larger uncertainties, but agree within the errors 
with the bcc results. r—■ 

For S = 1/2, the series 0(/3^^) of Ref.N yielded the estimates (without indication of error) a|“(l/2) = —0.11, 


^and = —0.13, which agrees closely with ours, 
obtained estimates of and for S = 1/2, 1 , 2, 


also quoted in Table tX and in 


which is 10 % larger than ourf 
Using the same series, Ref. I 
good agreement with ours. 

Our central estimate of the ratio a^/a^ is somewhat smaller than our prpz.ious estimateE^ a^/a^ = 0., 


shorter S' = 1/2 series, than the old HT estimates a^/a-^ = 0.83(5) of Ref.EII and a^/aj-f= 0.85 of Refs 


without indication of error), but it is larger than the HT estimate 0.71(7) of 
a^/a^ = 0.70(2), based on the fee series 0(/3^^) for general spin tabulated in Ref.| 



based on 
(reporter'' 


and the earlier estimateEj 
We should also mention the 


estimate a^/a^ = 0.65(5) obtained by RG in the peiiturbative fixed-dimension approach at sixth ordenHj. The 


e-expansion scheme (extending to second order) yieldedlHJ the estimate 
Fiaally, let us note that our results confirm the observations of Refs.HE 


F^tb65. 


and the arguments presented in 


Ref.Ed that the amplitudes of the leading CCS have a negative sign, both for the susceptibility and for the correlation 
length, in the case of the spin- 1/2 Ising model, on the sc and the bcc lattices. 


XIV. CONCLUSIONS 

For the Ising models of general spin S, on the sc and the bcc lattices, we have produced extended HT expansions 
of the nearest-neighbor correlation function, of the susceptibility, of the second correlation moment and of the second 
field derivative of the susceptibility. 

Our procedure of series analysis differs somewhat from the most traditional ones, but leads to completely consistent 
conclusions. At least for the models studied here, we are also confident that it yields very accurate direct estimates 
of the various critical parameters. Our updated results: 7 = 1.2371(1), i/ = 0.6299(2), 74 = 4.3647(20), gr = 1.404(3) 
and = 0.2668(5) are in good agreement with the latest calculations by other approximate methods, including 
the perturbative field theoretic RG approaches. At the same time, our new series data have proven to be sufficiently 
rich that we can obtain fairly tight checks of the conventional expectations about hyperscaling and universality, with 
regard both to the spin S and to the lattice structure. 
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TABLE I. Orders of high-temperature expansions, published (or obtainable from published data) before our work, for the 
nearest-neighbor correlation function G'(/3; S), for the susceptibility x{(^! *5')) for the second moment of the correlation function 


/r 2 (/ 3 ; S) and for the second field-derivative of the susceptibility S), in the case of the Ising models with general spin S. 


Observable Lattice; On 

te 

■[Ref.] Lattice; Ore 


Ref.] Lattice; Ore 


r[Ref.] 

x(/3; S) sc; 15 

fj ,2 {0; S) sc; 15 

Xi{(S-,S) sc; 14 

G((3-S) 

ZL 

in 

bcc; 21 
bcc; 21 
bcc; 10 

m 

in 

E 

fee; 14 
fee; 14 
fee; 13 
fee; 14 

U 

u 

u 



TABLE 11. Estimates of the critical inverse temperatures for the spin-5 Ising models on the sc and the bcc lattices, 
obtained in this work by the modified-ratio method biased with the leading correction-to-scaling exponent 9. The sensitivity 
of the estimates to the bias value of 9 is characterized by the derivatives of (3c- Eor comparison, we have also reported some 


results obtained by simulation methods or from the analysis of shorter series in the recent literature. 




S=l/2 

S=1 

S=3/2 

S=2 

S=5/2 

S=3 

S —CXD 

/3f(S) 

a/3*‘=(s)/ae 

0.221655(2) 

7■10"® 

0.312867(2) 

4•10“® 

0.368657(2) 

4■10“® 

0.406352(3) 

4■10“® 

0.433532(3) 

6■10“® 

0.454060(3) 

6■10“® 

0.601271(3) 

6■10“® 

/3c“(S) 

/d0 

0.1573725(10) 

3■10“® 

0.224656(1) 

10“® 

0.265641(1) 

-3 ■ 10“® 

0.293255(2) 

-6 ■ 10“® 

0.313130(2) 

-6 ■ 10“® 

0.328119(2) 

-6 ■ 10“® 

0.435085(3) 
-6 ■ 10“® 

/3r(s)E 


0.2216546(10) 







/3r(s)E 


0.221659555) 







/3r(s)b 

Jl 0.221655(1) 

WW) 

U 

0.157374(7) 

0.224657(4) 


0.293258(6) 



0.435089(11) 

/3r°(s)l 


0.157373 

0.224654 


0.293255 





TABLE III. Estimates of the critical exponents 7 and v for the spin-5 Ising model series on the sc and the bcc lattices 
obtained from extrapolation of the modified-ratio (MR) approximant sequences, (defined by eq.(|^)). We have also reported 
the estimates obtained in this paper by second-order inhomogeneous differential approximants (DA) biased with (3c- 


Exponent 

II 

S=1 

S=3/2 

S=2 

S=5/2 

S=3 

S—oo 

7 ““(S)(MR) 

1.2375(6) 

1.2378(7) 

1.2371(8) 

1.2367(10) 

1.2364(10) 

1.2363(10) 

1.2359(15) 

d'i‘'=(S)/d0 

0.016 

0.009 

0.007 

0.006 

0.006 

0.006 

0.006 

7 ''“(S)(MR) 

1.23742(20) 

1.23730(16) 

1.23710(3) 

1.23699(10) 

1.23694(10) 

1.23691(10) 

1.23685(15) 

a7'’“(s)/ae 

0.012 

0.003 

-0.0005 

-0.002 

-0.003 

-0.004 

-0.006 

7”VS')(DA) 

1.238(2) 

1.239(2) 

1.240(3) 

1.239(2) 

1.239(2) 

1.239(2) 

1.238(2) 

7 '’“(S)(DA) 

1.2378(8) 

1.2385(15) 

1.2370(4) 

1.2365(4) 

1.2366(4) 

1.2366(4) 

1.2367(4) 

!/““(S)(MR) 

0.6277(30) 

0.6279(30) 

0.6283(20) 

0.6285(20) 

0.6286(20) 

0.6286(20) 

0.6288(20) 

diy‘'‘{S)/d0 
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0.014 

0.010 

0.0095 
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0.008 

0.006 

i/'’“(S)(MR) 

0.6283(20) 

0.6294(8) 

0.6297(6) 

0.6298(6) 

0.6299(6) 

.6299(6) 

0.6301(5) 


0.006 

0.001 

-0.001 

-0.0015 

-0.002 

-0.0022 

-0.003 

!v““(S)(DA) 

0.632(2) 

0.631(1) 

0.631(1) 

0.631(1) 

0.631(1) 

0.631(1) 

0.631(1) 

i/'’“(S)(DA) 

0.6308(10) 

0.631(1) 

0.6299(4) 

0.6295(4) 

0.6296(4) 

.6295(4) 

0.6294(6) 


TABLE IV. Estimates of the exponents 7 and v obtained in the recent literature by various kinds of analyses of shorter 
high-temperature series, by MonteCarlo methods or by renormalization-group methods. The estimates marked with an asterisk 
are obtained by procedures implying or assuming universality. 


Exponent 


HT 


ft 


MonteCarlo 


RG e-exp. 


RG fixed-D exp. 


1.237(2) 

1.2385(25) 



1.2385(15)lfi 

1.2367(ll[lH 

□ 

1.2378(6)*|a 

7 

1.2395(4) ® 

1.2378(6)3 

1.2373(i; 

1.2371(4)ba 

1.237(2)*Sn 

1 . 237257 )^ 

1.2380(50)*0 

1.2396(13)’€3 


0.6300(15)* 

0.6305(15)[ 

0.632(1)* 

0.6311(3)f 

0.6300(2)1 

0.63002(23)1 


0.6296(7)1 

0.6301(8)*! 

0.6303(6)1 

0.6294(10)1. 


0.6305(25)*! 


0 


0.6301(5)*t 

0.6304(13)* 
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TABLE V. Estimates of the critical exponent 74 from the spin-S' Ising model high-temperature series on the sc and the bcc 
lattices obtained in this work by the modified-ratio (MR) method directly biased with the leading correction-to-scaling exponent 
9 following eq.(|3^) or by second-order differential approximants biased with fdc- For comparison, we have also reported some 
results obtained by simulation methods or from shorter series in the recent literature. 


Exponent S—1/2 S=1 S=3/2 S=2 S=5/2 S=3 S=oo 


7r(g) (DA) 

4.372(8) 

4.368(8) 

4.369(8) 

4.369(8) 

4.368(8) 

4.367(4) 

4.366(4) 

7r(g) (MR) 

dii'=(s)/ae 

4.3703(12) 

0.027 

4.3662(16) 

0.023 

4.3671(9) 

0.015 

4.3683(9) 

0.009 

4.3691(9) 

0.004 

4.3697(9) 

0.002 

4.3719(2) 

-0.007 

74 “°° (S) (DA) 

4.376(20) 

4.3666(10) 

4.3638(10) 

4.3629(10) 

4.3632(10) 

4.3631(10) 

4.3631(10) 


(MR) 
“'’“(St 


dll 


96 


4.3696(3) 

0.02 


4.3690(11) 4.3702(12) 4.3714(12) 

-0.008 -0.014 -0.023 


4.3722 4.3729(12) 4.3749(10) 

-0.028 -0.032 -0.043 


71° (g 

, bcc 


4.361(8) 


74 




4.366(6) 


74"°°(S)I 


4.370(14) 


TABLE VI. Estimates of the critical amplitudes A^{S) of the specific heat C^{P, S), C’^(5') of the susceptibility iS), 

f^{S) of the correlation length and C*{S) of xtiP'jlS) for the spin-S Ising models on the sc and the bcc lattices. 

They are obtained by differential approximants biased with the critical inverse temperatures reported in Table II and with 
the critical exponents estimated in this work. For comparison, we have also reported a few estimates obtained by simulation 


methods or from shorter series in the recent literature. For some of them, no indication of error is available. 




S=l/2 

S=1 

S=3/2 

S=2 

S=5/2 

S=3 

S = CXD 

A“(S) 


04(11 

1.80(3) 

2.00(2) 

2.09(1) 

2.15(2) 

2.18(2) 

2.28(1) 

aA“(5)/d/3; (g) 

640. 

980. 

170. 

550. 

300. 

750. 

440. 

9A“(S)/aa 

-20. 

-26. 

-28. 

-32. 

-32. 

-32. 

-33. 

A'’“(S) 


1.302(6) 

1.732(6) 

1.911(6) 

2 . 000 ( 6 ) 

2.051(6) 

2.082(6) 

2.171(6) 

9A““(S)/a/3f=(S) 

800. 

1000. 

950. 

900. 
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870. 

830. 

3A“'=(£ 

Uda 

-19. 

-26. 

-29. 

-30. 

-31. 

-32. 

-33. 

A‘“={S)b 


1.464(7) 







A“(S)Ei 

Li - 

1.368(7) 







A'’“(sl 

2U 

1.431(9) 







C“(S) 


1.127(3) 

0.682(1) 

0 . 545 ( 1 ) 

0.482(1) 

0.443(1) 

0.4184(6) 

0.3073(2) 

9C“(g)/a/3“(g) 

1900. 

1800. 
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600. 

220. 

150. 

dC“‘^{S)/di 

-9. 

-8. 

-7. 
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-5.5 

-5. 

-3. 

C'’“(S) 


1.042(1) 

0.622(1) 

0.4967(3) 

0.4379(3) 

0.4045(4) 

0.3826(4) 

0.2817(4) 

9C“°(S)/a/3c°°(g) 

3900. 

-1300. 

170. 

750. 

500. 

480. 

310. 



- 9 . 

-5. 

-4.5 

-3.6 

-3.1 

-3. 

-2.2 

C“(S)1: 

2L 

1.1025(10) 







C“(S)E 

i2 

1.093(13) 







C'’“(S) 


1.026 

0.620 


0.4346 




C“°(S)[ 


1.0312(10) 







/“(g) 


0.506(1) 

0.458(1) 

0.443(1) 

0.436(1) 

0.432(1) 

0.430(1) 

0.423(1) 

a/“(g)/a/3r(g) 

290. 

280. 

200. 

180. 

400. 

120. 

120. 

df‘^{s)/du 

-4. 

-4. 

-3.5 

-4. 

-4. 

-3. 

-3.5 

/““(g) 


0.4686(4) 

0.4262(8) 

0.4112(4) 

0.4047(4) 

0.4013(4) 

0.3992(4) 

0.3937(4) 

9/"“(g)/a/3““(g) 

500. 
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-5. 

- 3 . 
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-3. 
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-3. 

/“(g) 

1 

0.496(4) 







/“(g)“ 

‘ 

0.501(2) 







/“(g)ffi 


0.5192 







/““(g:E 


0.4590(1) 







/““(g)L 

I 

0.46821 

0.42605 


0.4038 




Cl‘^(S) 


3.87(1) 

1.05(1) 

0.606(1) 

0.456(2) 

0.375(1) 

0.327(1) 

0.169(1) 

gC“(S)/a/3“(g) 

5000. 

3200. 

2900. 

2700. 

800. 

-1300. 

350. 

dcms)/dii 

-20. 

-9. 

- 7 . 

-6. 

-5. 

-1.5 

-1.5 

Gl'=‘(S) 


3.410(8) 

0.912(3) 

0.523(1) 

0.3884(6) 

0.3230(5) 

0.2847(5) 

0.1478(5) 

dCl‘^‘=(S)/d/9l‘^‘=(S) 

15000. 

8000. 

460. 

630. 

470. 

450. 

190. 

ocr°(|S; 

\/dlA 

-24. 

-9.5 

-3. 

-2.5 

-2. 

-2. 

-1. 

cr(g)= 

2 

3.70(3) 








C|°°(S)Liil 3.236(2) 
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TABLE VII. Estimates of the universal renormalized coupling Qr using the auxiliary function z^{(3\S) of eq.(^^ for the 
spin-S Ising model series on the sc and the bcc lattices. They are obtained by differential approximants biased with the 
modified-ratio estimates of the critical inverse temperatures reported in Table II. For comparison, we have also reported other 
estimates obtained by simulation methods or from shorter series in the recent literature. (For some of them, no indication 
of error is available.) The values marked with an asterisk have been obtained either by renormalization-group methods or 
by high-temperature methods which assume universality and therefore they refer to the Ising universality class although, for 
simplicity, they are reported in the column of the 5 = 1/2 results. 


S=l/2 S=1 S=3/2 S=2 S=5/2 S=3 S=oo 

g“‘^{S) (DA) 1.40(1) 1.410(6) 1.404(6) 1.414(10) 1.415(10) 1.414(10) 1.412(10) 



(DA) 1.409(4) 1.404(3) 1.40l[/*( l-400[/5( l-400}/®' 


bid 1.401(8) 

art 

1.403(7) 


3 1.459(9) 

9T-\h 

1.402(2)* 

gr- 

1.40* 

3.0 

1.411(4)* 


TABLE VIII. Estimates of the universal quantity using the auxiliary function F of eq.(^3|) for the spin-S Ising model 
series on the sc and the bcc lattices. They are obtained by differential approximants biased with the critical inverse temper¬ 
atures reported in Table II and with the value of v obtained in this study. We have also reported some estimates obtained 
by other methods or from shorter series in the recent literature. The estimates marked with an asterisk are obtained by 
renormalization-group methods or by high-temperature methods that assume universality and therefore they refer to the Ising 
universality class although, for simplicity, they are reported in the column of the 5=1/2 results. 


S=l/2 


S=1 


S=3/2 


S=2 


S=5/2 


S=3 


S=oo 



TABLE IX. Estimates of the universal ratio a^/a^ for the spin-5 Ising model series on the bcc lattice. For comparison, we 
have also reported some estimates obtained by renormalization-group methods. (For some of them, no indication of error is 
available.) Although they refer to the Ising universality class, for simplicity, they are reported in the column of the 5=1/2 


results. 





S=l/2 

S=1 

S=3/2 

S=2 

S=5/2 

S=3 

S—oo 

at“(S) 


-0.129(3) 

-0.0363(10) 

0.0079(8) 

0.0307(10) 

0.0436(10) 

0.0515(10) 

0.0742(20) 

4“(S) 


-0.100(4) 

-0.0279(15) 

0.0061(6) 

0.0233(10) 

0.0331(20) 

0.0390(20) 

0.0560(20) 

a^“(S) 


-0.119 

-0.034 


0.023 




4“(S) 


-0.1085 

-0.033 


0.0225 




a^“(S) 


-0.13 

a^‘=‘=(S) 

21 

-0.11 



i““(S) 0.78(5) 

0.77(6) 

0.77(16) 

0.76(6) 

0.76(6) 

0.76(5) 

0.76(5) 

a£/oj 

'iik . 

0.85 



0.65(5)* 


VIA 

0.65* 
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FIG. 1. The “normalized” modified-ratio approximant sequences (/3*’““(5'))„/A'^*’““(S') for the critical inverse temperature of 
the spin S Ising models on the bcc lattice, plotted versus l/n^+®, with 6 = 0”'^^ = 0.504. They are obtained from eq.(|2^) 
using the coefficients c^'^{S) of the susceptibility series for the bcc lattice. In order to fit into a single figure the sequences 
for different values of the spin, each sequence has been normalized by the average of the critical inverse temperatures 

obtained extrapolating separately the even and the odd subsequences. We have indicated by continuous lines the extrapolants 
of the odd subsequences, based on the last odd pair of approximants, while the dashed lines indicate the extrapolants of the 
even subsequences, based on the last even pair of approximants. 


FIG. 2. The same i 
from the coefficients cii 



(S) of the susceptiblity series for the sc lattice. 


"(5) formed 
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1/n® 

FIG. 3. The modified-ratio approximant sequences of the susceptibility critical exponents for various values of 

the spin S, plotted vs. 1/n®, with 0 = = 0.504. They are obtained using eq.(|2^) from the susceptibility series coefficients 

c^'^{S). For each value of the spin, we have indicated by a continuous line the extrapolation to large n of the sequence, linearly 
in 1/n®, based on the last odd pair of terms {(7*’“(5'))23, ( 7 *’'^'^(S)) 25 }. A dashed line indicates the extrapolation based on the 
previous odd pair {( 7 *’‘'‘'(S')) 2 i, (7'’“(5'))23} 



FIG. 4. The sequences of the extrapolations of the successive odd pairs of modified-ratio approximants (7*’'^'^(S))n vs. 
l/n^+®.The continuous lines are only guides to the eye. 
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FIG. 5. Same as in Fig. ^ but for the modified-ratio approximant sequences of the correlation-length critical 

exponent, as obtained from the expansion of S). 



0.0 0.002 0.004 0.006 0.008 0.010 0.012 0.014 

; ( 1 + 3 ) 

1—1 1 / n 

FIG. 6. Same as in Fig. m but for the sequences of the extrapolations of the successive odd pairs of modified-ratio 
approximants vs. 
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FIG. 7. Same as in Fig. but for the modified-ratio approximant sequences ( 7 ®'^(S'))n of the susceptibility critical exponent 
as obtained from S) using eq.(Gffl. 



FIG. 8. Same as in Fig. but for the modified-ratio approximant sequences {L'“^{S))n of the correlation-length critical 
exponent, as obtained from 
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FIG. 9. The directly biased modified-ratio approximant sequences (7 *’“(S'))to plotted vs. They are obtained from 

eq.(^), using as a bias the value of 6 in order to reduce the influence of the confluent corrections to scaling. The continuous 
lines are only guides to the eye. 


o 

eg 



FIG. 10. Same as Fig. but for the directly biased modified-ratio approximant sequences plotted vs. 
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0.0 0.0002 0.0004 0.000<5 0.0008 0.0010 

FIG. 11. The directly biased modified-ratio approximant sequences (j 4 ‘^‘^(S))n for the critical exponent of S) plotted 

vs. l/n^"*"®. In order to keep the figure readable we have indicated only the extrapolations of the odd approximant subsequences. 



0.0 0.0002 0.0004 O.OOOd 0.0008 0.0010 

. ( 2 + 8 ) 

|—I 1/n 

FIG. 12. Same as in Fig. 1^, but for the directly biased modified-ratio approximant sequences (74'^(5'))n for the critical 
exponent of xTiP'i S) plotted vs. l/n^+®. 
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FIG. 13. Approximant sequences S 2 ))n for the correction-to-scaling exponent 9 as obtained using eq.(^5|) for fixed 

5i = 1/2 and S 2 = 5/2,3,7/2,4,5,00. The symbols refer to the values of S 2 . The sequences are plotted vs. l/n^+®, with 
9 = 0.504. In order to keep the figure readable we have indicated only the extrapolations of the odd approximant subsequences. 



a 

X 

FIG. 14. Highest order simplified-differential approximants of the effective exponent 7e/‘/(/3; S) of the susceptibility S) 

as defined by eq.(^^. For each value of the spin S the effective exponent is plotted vs. = (1 — fi/{S)Y. 

As indicated by the symbols attached to them, the curves refer, from the highest downwards, to the spin values 
S = 1/2,1, 3/2, 2, 5/2, 3, 7/2, 4, 5, 00 . 
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FIG. 15. Highest order simplified-differential approximants of the effective exponent v^Yf{l3\S) of the correlation length 
S) as defined by eg.. For each value of the spin S the effective exponent is plotted vs. = (1 — f3/(S))^. 

As indicated by the symbols attached to them, the curves refer, from the highest downwards, to the spin values 
S = 1/2,1, 3/2, 2, 5/2, 3, 7/2, 4, 5, oo. 


o 



fl 

X 

FIG. 16. Highest order simplified-differential approximants of the effective exponent S) of the susceptibilty S) 

as defined by eq.(^). For each value of the spin S the effective exponent is plotted vs. = (1 — (5/131’^{S)Y. 

As indicated by the symbols attached to them, the curves refer, from the highest downwards, to the spin values 
S = 1/2,1, 3/2, 2, 5/2, 3, 7/2, 4, 5, oo. 
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FIG. 17. Highest order simplified-differential approximants of the effective exponent of the correlation length 

as defined by eg. (|3^) . For each value of the spin S the effective exponent is plotted vs. = (1 — P/(3"'^{S))^. 

As indicated by the symbols attached to them, the curves refer, from the highest downwards, to the spin values 
S = 1/2,1, 3/2, 2, 5/2, 3, 7/2, 4, 5, oo. 


CO 
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T 

FIG. 18. Highest order simplified-differential approximants of the effective exponent 74e//(/3; S) computed from S)- 

For each value of the spin S the effective exponent is plotted versus the corresponding reduced inverse temperature 
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CO 



8 

X 

FIG. 19. Same as Fig. |^, but for the effective exponent 'yllffiP', S) computed from S). 



FIG. 20. Highest order simplified differential approximants of the effective dimensionless renormalized coupling constant 
gr{f3',S), as obtained from the bcc lattice series. The effective coupling is computed from the auxiliary function S') 

defined in eq.(^ and is plotted vs. = 1 — /l//3c“(S). 
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FIG. 21. Highest order simplified differential approximants of the effective dimensionless renormalized conpling constant 
grr(/3; S) as obtained from the sc lattice series. The effective coupling is computed from the auxiliary function S) defined 

by eq.(^^ and is plotted vs. = 1 — P/l3c‘^{S). 
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